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Abstract 

We present the first calculation of the one-loop QCD corrections for the decay of 
an off-shell vector boson with vector couplings into two pairs of quarks of equal or 
unequal flavours keeping all orders in the number of colours. These matrix elements 
are relevant for the calculation of the next-to-leading order 0{o? s ) corrections to four 
jet production in electron-positron annihilation, the production of a gauge boson ac- 
companied by two jets in hadron-hadron collisions and three jet production in deep 
inelastic scattering. We compute the interference of one-loop and tree level Feynman 
diagrams, and organise the matrix elements in terms of combinations of scalar loop in- 
tegrals that are well behaved in the limit of vanishing Gram determinants. The results 
are therefore numerically stable and ready to be implemented in next-to-leading order 
Monte Carlo calculations of the e + e~ — > 4 jet, e ± p — > + 3 jet and pp — > V + 2 jet 
processes. 



1 Address after 1 January 1997, Rutherford Appleton Laboratory, Chilton, Didcot, Oxon, 0X11 OQX, 
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An important ingredient in the confrontation of hadronic jet data from high energy par- 
ticle collisions with theoretical calculations in perturbative QCD is the calculation of parton 
level cross sections beyond leading order. This is so for a variety of reasons. First, the 
dependence on the unphysical renormalisation and factorisation scales is usually reduced. 
Second, the addition of more particles into the final state makes the calculation more sen- 
sitive to the details of the experimental jet algorithm. Third, the final state phase space is 
enlarged so that lowest order constraints on the jet transverse energy or rapidity are relaxed. 
Together, these improvements make for a more stable calculation that is better matched to 
the experimental data. Finally, the presence of large infrared logarithms can be detected 
and the appropriate resummations performed. In this paper, we present the first calculation 
of the one-loop matrix elements for a virtual gauge boson decay into two pairs of equal or 
unequal flavour quarks keeping all terms in the number of colours. These matrix elements 
are one of the ingredients of the next-to-leading order calculation relevant for the processes 

• e + e~ — > 4 jets 

• e ± p — >■ e ± + 3 jets 

• pp — * V + 2 jets with V = W or Z, 



which have been observed in abundance in experiments at the major accelerators, LEP at 
CERN, HERA at DESY and the TEVATRON at Fermilab. In electron-positron collisions, 
applications include more precise measurements of the colour casimirs Ca/Cf and T F /Cp 
and a reduction of the theoretical error in the extraction of the WW cross section near 
threshold. 

In recent years there have been several significant technical advances in one-loop Feynman 
diagram calculations, some relying on insight from supersymmetry and others inspired more 
directly by the superstring Most important has been the continuation of the one- loop 
scalar pentagon integral in 4-dimensions [F2] to D = 4 — 2e dimensions J3], ^, At the 
same time, methods utilised in tree level matrix element calculations, colour decompositions 
and helicity have been taken over. To make best use of the helicity method, relations have 
been provided between conventional dimensional regularisation, where all polarisations and 
spinors are treated in D- dimensions, and other schemes where parts of the calculations are 
performed in 4-dimensions || [7]]. With these technological improvements, compact one- 
loop helicity amplitudes for processes involving five partons have been presented by Bern, 
Kosower and Dixon [| and Kunszt, Signer and Trocsanyi JTU], [TT], |T2"|| , while recently, 
one of the helicity amplitudes for the leading colour part of the Z — > qqgg process has 
been evaluated ||13|| . Recently, Dixon and Signer have reported first numerical results for 



the leading colour contribution to the e + e — > 4 jet rate at next-to-leading order fL"4"f . As 



expected, the agreement between theory and experiment is significantly improved. Our 



1 



results presented here can therefore provide a numerical check of the one-loop leading colour 
contribution from the four quark process.^ 

For our calculation involving a massive external gauge boson, we choose to compute the 
interference between the tree-level and one-loop amplitudes directly. The resulting 'squared' 
matrix elements involve only dot products of the particle four-momenta in addition to the 
normal logarithms and dilogarithms encountered in one-loop Feynman diagram calculations. 
As a consequence, we can immediately eliminate all integrals more complicated than a scalar 
pentagon, although some box tensor-integrals with three loop momenta remain. We have 
kept all quantities in .D-dimensions as in conventional dimensional regularisation, but have 
checked that the 't Hooft-Veltman scheme [IB] gives identical results |], 0]. Throughout we 



reorganise the tensor loop integrals as combinations of scalar integrals using a momentum 
decomposition jl7| (rather than the reciprocal momentum basis of [[18], || |12|). The usual 
problem with this approach is the occurence and proliferation of Gram determinants (A) in 
the denominator which both increases the size of the algebraic expressions and introduces 
numerical instabilities as A — > 0. We have avoided this by introducing combinations of 
scalar loop integrals in groups that are finite in the limit that the Gram determinants vanish 



H [10], [12], [H|. The matrix elements, written in terms of these finite functions, are then 
obviously well behaved in the same limits. These groupings are derived by differentiating 
the scalar integrals with respect to the external kinematic factors or by considering them 



in higher dimensions [H|. They combine the dilogarithms, logarithms and constants from 



different scalar integrals in an extremely non-trivial way and, in doing so, reduce the size of 
the expressions considerably. However, these finite combinations are not linearly independent 
(unlike the raw scalar integrals) and by suitable rearrangement of the polynomial coefficients, 
one function may be transformed into another, and some ambiguity in the presentation of 
the final answer remains^. To generate and simplify our results, we have made repeated use 
of the algebraic manipulation packages FORM and MAPLE. 

We choose to work with all particles in the physical channel corresponding to e + e~ 
annihilation where all quarks are in the final state, 7* — > qqQQ. The momenta are labeled 

as, 

7*(Pi234) -> q(pi) + q{p2) + Q(p 3 ) + Q(P4). (1) 

Using momentum conservation, we can systematically eliminate the photon momentum in 
favour of the four massless quark momenta. 

For quarks q and Q, the electric charges are denoted e q and tQ while the colours of the 



2 Since this paper was completed, Bern et al. |L5| have presented analytic results for all of the helicity 
amplitudes for the V — * QQlQ process. Their approach is completely different from ours and, using the 
known relations between amplitudes calculated in dimensional regularisation and the 4-D helicity scheme 
[^] a numerical comparison will provide a powerful check of our results. 

^3ven when only scalar integrals are present, the form of the polynomial coefficient depends on how the 
Gram determinants have been cancelled. 
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quarks are labeled by q, i — 1, . . . , 4. The colour structure of the matrix element at tree-level 
(n = 0) and one-loop (n = 1) is rather simple and we have, 



" 2 UnJ 



KcA 3 c 2 U B) (1,2) +4 n) (3,4) + % (4 n) (l,4) + 4 n) (3,2; 



- Wcc (4 n) (l, 2) + 4 n) (3, 4)) + 5 qQ (4 n) (l, 4) + 4 n) (3, 2))) j , (2) 



AT 

where N is the number of colours and 8 q g = 1 for identical quarks and zero otherwise. The 
arguments of Ai indicate which quark line is attached to the vector boson and hence which 
quark charge that function is proportional to. At lowest order, 

A?\i,j) = A$ ) (i,j), (3) 

while at one-loop we find, 

A?(i,j) = NJ$\iJ)-±(2A%\iJ)+A$(i,j)), (4) 
= N(A^(iJ)-A^(iJ))-^(A ( i\iJ)+A^(i,j)). (5) 



The functions A)p(i,j), a = A, B,C represent the contributions of the three gauge invariant 
sets of Feynman diagrams shown in Fig. 1 where the photon couples to the quark-antiquark 
pair i, j. The set of diagrams containing closed fermion triangles do not contribute by Furry's 
theorem. Note that we choose to include the contribution from the fermion loop in the fourth 
diagram of Fig. 1(c) which is proportional to the number of flavours Np in the leading colour 
part Ac- 



The squared matrix elements are straightforwardly obtained at leading order p0[ . Be- 
cause the matrix element contracted with the photon momentum vanishes by current conser- 
vation, the sum over the polarisations of the gauge boson may be performed in the Feynman 
gauge, 

E ^ = ~9" U - (6) 

spins 

Hence, 

2^4 



£ |^(0)|2 = iA(iV2_l) 



spins 

X 



4 

T(l, 2; 1, 2) + T(l, 2; 3, 4)) + 6 -f (V(l, 2; 1, 4) + T(l, 2; 3, 2)) j 

1 «_> 3; 2 <-> 4) + 5 qQ {2 <-> 4) + 5 qQ (l <-> 3), (7) 
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(a) 






Figure 1: The classes of Feynman diagrams relevant for the different colour structures. The 
solid circle indicates the possible positions for attaching the off-shell photon to the quark- 
antiquark pair Group (a) contributes to A^\i,j), (b) to Ag\i,j) and group (c) to the 
leading colour amplitude A^\i,j). Diagrams with self-energy corrections on the external 
lines are zero in dimensional regularisation and have been omitted. 
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where, 



T(i,j;k,l)= £ \A<M(i,j)A®(k,l)\. 

spins 



(8) 



The relevant 'squared' matrix elements are the interference between the tree- level and 
one-loop amplitudes, 



J2 2\M m M {1 



spins 
X 



£ c (l,2;l,2) + £ c (l,2;3,4) 
1 



N 2 



2C A (1, 2; 1, 2) + 2C A (1, 2; 3, 4) + C B (1, 2; 1, 2) + £ B (1, 2; 3, 4) 



+5, 



qQ 



N 



C c (l, 2; 1, 4) + £ c (l, 2; 3, 2) - £ A (1, 2; 1, 4) - £ A (1, 2; 3, 2) 



±( 

N 3 V 



£ A (1,2;1,4) + £ A (1,2;3,2) + £ B (1,2;1,4) + £ B (1,2;3,2; 



+ (1 ^ 3, 2 ^ 4) + 5, Q (2 ^ 4) + 5, Q (1 <-> 3), 



with, 



(9) 
(10) 

(11) 



spins 

Using the symmetry properties of the Feynman diagrams, we find that, 

C a (l, 2; 3, 2) = £ Q (1, 2; 1, 4)( Pl ^ p 2 ,p 3 p 4 ), 

so that the 'squared' matrix elements are described by 9 independent £ a . 

In computing the 'square', the loop momentum always appears as either C 2 or E.p. We 
therefore systematically reduce the number of loop momenta appearing in the numerator by 
simple rewriting, 



■P 



(£ + P f -i 2 -p 2 



£ 2 (£ + p) 2 2£ 2 (£ + p) 2 

The only time that we are unable to do this occurs when the loop momentum is contracted 
with the 'wrong' momentum, i.e., one that does not appear in any of the propagators. For the 
pentagon diagrams, all momenta are involved in the propagators and we can always rewrite 
the integral as a combination of box tensor-integrals. Dividing through in this way, we find 
box (triangle) tensor-integrals with at most three (two) loop momenta in the numerator. 

The simplified tensor integrals can now be expressed in terms of scalar integrals using 
a momentum decomposition. We choose the natural momentum as in [17] rather than 
the reciprocal momentum basis of JIR |T^]. The coefficients of the scalar integrals are 



obtained by solving systems of linear equations and the usual problem is the occurrence 
of Gram determinants. Typically, for a rank n tensor m-point diagram, we can obtain an 
(m — 1) x (m — 1) Gram determinant raised to the nth power in the denominator. Since the 
physical process does not usually possess singularities near the edge of phase space where 
the Gram determinants vanish, one obtains large cancellations between terms and potential 
numerical inaccuracies. For 'squared' one- loop matrix elements the number of terms is 
expected to be large, thereby exacerbating the numerical problem. However, since these 
Gram determinant singularities are unphysical, it should prove possible to arrange the scalar 
integrals in combinations that are finite in the limit that the Gram determinant vanishes. 
The matrix elements, written in terms of these finite combinations, are then obviously finite 
and numerically stable in the same limits. The explicit results of Bern, Kosower and Dixon 
|| ^] and Kunszt, Signer and Trocsanyi [10, 12] for the five-parton one- loop amplitudes 



already show the introduction of such functions. For example, Bern, Kosower and Dixon 
introduce the functions 

1 — r 1 — r 

which are finite as r — > 1. Furthermore, the helicity coefficient of these functions is usually 
rather simple so one might expect that the 'squared' matrix element coefficient is not enor- 
mous. We have therefore assembled a collection of finite functions relevant to the process 
at hand jnj and organised the C a as a linear combination of these functions where the 
coefficients are polynomials of the generalised Mandelstam invariants, 

Sij = (Pi + Pj) 2 , s ijk = (pi + pj + p k ) 2 , s m = (pi + pj +p k + Pi) 2 - 

As mentioned earlier, these fuctions are natural in the sense that they are obtained by 
differentiating the scalar integral with respect to the external kinematic parameters or by 
evaluating the integral in higher dimensions. However, these finite combinations are not 
linearly independent (unlike the raw scalar integrals), and by explicitly cancelling part of 
the polynomial coefficient against the Gram determinant it is possible to shift terms from 
one function to another. Therefore an ambiguity in the presentation of the final answer 
remains. 

As mentioned earlier, we use dimensional regularisation and work in D = 4 — 2e dimen- 
sions. It is straightforward to remove the infrared and ultraviolet singularities poles from 
the C a since they are proportional to the tree-level amplitudes, 

£a(1,2;*,j) = (+^-^-^ + ^)r(l,2;i,i)+£ A (l,2;i,j), (12) 
C B (WJ) = (-^-^-^)r(l,2;z,j) + ^(l,2;z,j), (13) 
£ c (l,2;i,j) = [-— — + -—-—— —jr{l,2]i,j)+C c {l,2]i,j), (14) 
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where we have introduced the notation, 



v " = {—,) r ( i-2e) • (15) 

In physical cross sections, this pole structure must cancel with the infrared poles from 
the 7* — > qqQQ + g process and those generated by ultraviolet renormalisation. This pole 
structure is in agreement with the expectations of ref. ETJ and reproduces that given in 



11, 12 



The finite functions C can be written in the following symbolic way, 

C = J2Pi(s)U, (16) 

i 

where Pi{s) is a polynomial of the invariant masses s^- and the Lj are linear combinations of 
scalar loop integrals. Rather than have the Gram determinants present in the polynomial 



coefficients as in the conventional method |17], I8|, [| 0), they are absorbed into the L 



which are constructed to be well behaved in the limit that the Gram determinant vanishes. 
The advantage is better numerical stability while the disadvantage is that there are rather 
more functions (than the raw scalar integrals in 4 — 2e dimensions). Furthermore, since the 
functions are not linearly independent the polynomial coefficients are not unique. Typically, 
the Pi{s) contains O(30) terms (roughly the same size as the tree level result) while the 
summation runs over O(30) functions. Although the expreassions for the individual C a are 
in closed form, they are still rather lengthy and each contains several hundred terms. We 
have therefore constructed a FORTRAN subroutine detailing the P»(s) and evaluating the 
finite one-loop contribution for a given phase space point. This will act as an ingredient in 
calculating the full next-to-leading order e + e~ — > 4 jet rate once the necessary sub-leading 
colour contributions to the 7* — ► qqgg matrix elements are known via numerical combination 



with the tree level 7* — > 5 parton processes [ 22J, ^3], pM |. 



To illustrate the finite groupings Lj and the potential numerical instabilities, we show 
the functions associated with the triangle graph with three massive external legs. The scalar 
triangle integral with momenta P12 and P34 flowing out (and P1234 flowing in) can be written 



25. 



Lc (pi2,P34) = -i= Hog(a + a-) log (j^) + 2Li 2( a+ ) ~ 2Li 2 (a") j , (17) 



where Li 2 is the usual dilogarithm function and a ± are two roots of a quadratic equation 



a = , (18) 

^■§1234 
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and, 

A 3 = — S 1234 — S 12 — S u + 2S1234S12 + 2S1234S34 + 2s 12 S34. (19) 

The functions, 

Lcis(Pl2,P34) = ^-f2Si234Sl2S 3 4LCo(j5l2,P34) (20) 



Su(Sl23A + S34 - S U ) log ( ) + S 34 (Si234 + «12 ~ S 3A ) log 

V S12 / 



•§1234 



12 / V S34 

Lc 2 s(Pl2,P34) = ^-^2Si234Sl2S34LCi,s(Pl2,P34) - ^1234^12534 (21) 

1 / 2 / \ 1 / S 1234\ . 2 / M / s 1234 

-77 S 12 (Si234 + S34 - S12) log + S 34 (Si234 + S 12 ~ S34) log 

6 V V S12 / V S34 

also appear and are related to the scalar triangle integral in D = 6 — 2e and £) = 8 — 2e 
dimensions [19|. Throughout, the arguments of the functions describe the momenta flowing 
out of all but one of the vertices, while the final momentum is determined by momentum 
conservation. For triangle integrals it is also convenient to introduce the additional functions 
corresponding to integrals in 4 — 2e dimensions with Feynman parameters in the numerator 
p9|, for example, 



LciG»12,£>34) = ^-fsi 2 (Si234 - Si2 + S 3 4)LC (P12,P34) 

-(S1234 - Su - S 34 ) log (~f^j ~ 2S12 log (~^) J ' ( 22 ) 
Lc 2 (pi2,P34) = ^ 2s 34(Sl234 + «12 ~ S 34 )Lci (p 12 , p U ) + Sl 2 (Sl234 + S34 ~ S 12 )LCi (p 34 , p 12 ) 

-Si2S 34 Lc (pi2, P34) - S34 log (-J^j + S 12 + S34 - Sl234^j • (23) 

As A 3 — > 0, all of these functions are finite and combine dilogarithms and logarithms in a 
highly non-trivial way. For example, 

lim Lcis(pi 2 ,P34) = 77; " ( s? 2 ( s i234 + S34 - su) log (^^) 

A 3 ^0 12Si 2 S 34 Si234 V V S12 / 

+ S3 4 (Sl234 + S 12 - S 34) log {jj^j + 2 -5l2S 3 4Sl234^ • (24) 

By keeping terms proportional to A 3 (and higher), the limit can be approached with arbitrary 
precision. 
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Figure 2: The finite functions for the triply massive triangle graph with S1234 = 1 and 
S12 = 0.2 as a function of A 3 /A™ ax where A™ ax = — (S1234 — s^) 2 - The dashed lines show 
the approximate form for the function in the limit A 3 — > 0, retaining only the first term of 
the Taylor expansion. 



Numerically, the instability due to the Gram determinant is easy to see if we restrict 
ourselves to a specific phase space point, S1234 = 1, S12 = 0.2 and letting s 34 vary in such a 
way that A 3 — > 0. This corresponds to S34 — > 0.135 which is neither the soft nor collinear limit 
and there should be no kinematic singularity. The behaviour of the various triangle functions 
is illustrated in Fig. ^| as a function of As/A™* where A™ x = — (si23 4 — S12) 2 For a given 
numerical precision, acc, the numerical problems typically occur when A3/ A™ x ~ (acc) 1 ^ 
where N is the number of Gram determinants in the denominator of the function. Since the 
scalar integral already contains A3 1 ^ 2 , for an intrinsic numerical precision of roughly 10~ 14 , 
the Lei and Lcis functions break down at A 3 /A™ x ~ few x 10~ 10 while the LC2 and LC25 
functions with effectively 2.5 inverse powers of A3 break down at A^/A^ x ~ few x 10~ 6 . 
Other phase space points yield a similar behaviour. With adaptive Monte Carlo methods 



such as VEGAS [26], once the singular region is found, more and more Monte Carlo points 



are generated there and the result will be unreliable. However, for A 3 < 10 _5 A 3 nax we 
see that the approximate form obtained by making a Taylor expansion about A 3 = and 
keeping only the constant term works well. 

Similar groupings for the box and pentagon graphs (with their corresponding Gram 
determinants) can be obtained by considering the scalar integrals in higher dimensions or 



with Feynman parameters in the numerator This approach is easily generalisable to 

other processes with more general kinematics. 
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To summarize, we have performed the first calculation of the one-loop 'squared' matrix 
elements for the 7* — > qqQQ process keeping all orders in the number of colours. Throughout 
conventional dimensional regularisation has been employed as well as grouping the Feynman 
diagrams according to the colour structure. However, in order to deal with the Gram deter- 
minants, we have expressed the formfactors in terms of functions (Lj) which group together 
dilogarithms, logarithms and constants from the basic scalar integrals in a way that is well 
behaved as the Gram determinants vanish. These groupings are derived by differentiating 
the scalar integrals with respect to the external kinematic factors or by considering them in 
higher dimensions. The resulting expressions are still rather long, due partly to the number 
of functions and partly to the size of the tree-level matrix elements, but numerically sta- 
ble and ready to be implemented in next-to-leading order Monte Carlo calculations of the 
e + e~ — > 4 jet, e ± p — > e ± + 3 jet and pp — > V + 2 jet processes. 
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